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Abstract 

The semiclassical limit of the algebraic quantum inverse scattering method 
is used to solve the theory of the Gaudin model. Via Off-shell Bethe ansatz 
equations of an integrable representation of the graded osp{\\2) vertex model 
we find the spectrum of the N — 1 independents Hamiltonians of Gaudin. In- 
tegral representations of the iV-point correlators are presented as solutions of 
the Knizhnik-Zamolodchikov equation. These results are extended for highest 
representations of the osp(l\2) Gaudin algebra. 
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1 Introduction 



In integrable models of statistical mechanics |]T|, an important object is the 7£-matrix 
where A is the spectral parameter. It acts on the tensor product V 1 £g> V 2 for a 
given vector space V and it is the solution of the Yang-Baxter (YB) equation 

K 12 (\)K 13 (\ + fi)-Rw(ji) = n 23 (fi)TZ 13 (X + fM)K 12 (\) (1.1) 

in V 1 ® V 2 ® V 3 , where 7Zi 2 = 7Z <S>X, 7Z 23 = X® 7?., etc. and X is the identity matrix. 
If 7?. depends on a Planck-type parameter 77 so that 1Z(\,r)) = 1 + 77 r(A) + o(?y 2 ), as 
?7 — ^ 0, then the "classical r-matrix" obeys the classical YB equation 

MA), r 13 (A + /i) + r 23 (M)] + MA + m), MaO] = (1.2) 

Nondegenerate solutions of ( |1.2| ) in the tensor product of two copies of simple Lie 
algebra g , r\,(A) G g» <8> gj , i,j = 1, 2, 3, were classified by Belavin and Drinfeld [T2J. 

The classical YB equation interplays with conformal field theory in the following 
way: In the skew-symmetric case rji(— A) +ry(A) = 0, it is the compatibility condition 
for the system of linear differential equations 

d 

K—ty(z u ...,z N ) = Y, r ij( z i - z j)^( z ii ■■■iZn) (1-3) 

in Y complex variables Zi, Zn for vector- valued functions ^ with values in the tensor 
space = V 1 <E> • • • <8) and k is a coupling constant. 

In the rational case ||, very simple skew-symmetric solutions are known: r(A) = 
C 2 /A, where C 2 G g ® g is a symmetric invariant tensor of a finite dimensional Lie 
algebra g acting on a representation space V. The corresponding system of linear 
differential equations is the Knizhnik-Zamolodchikov (KZ) equation for the conformal 
blocks of the Wess-Zumino-Novikov-Witten (WZNW) model of conformal theory on the 
sphere [[|. Therefore, one can consider a KZ equation for each solution of the classical 
YB equation. 

The algebraic Bethe ansatz Q is the powerful method in the analysis of integrable 
models. Besides describing the spectra of quantum integrable systems, the Bethe ansatz 
also is used to find reasonably efficient expressions for the correlators ||. Various 
representations of correlators in these models were found by Korepin 0, using this 
method. 

Recently, Babujian and Flume [0 developed a method which reveals a link to 
the algebraic Bethe ansatz for the theory of the Gaudin model. Their method was 
confirmed by another different approach presented in the work of Feigin, Frenkel and 
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Reshetikhin ||. In the Babujian-Flume method the wave vectors of the Bethe ansatz 
equation for inhomogeneous lattice model render in the semiclassical limit solutions of 
the KZ equation for the case of simple Lie algebras of higher rank. More precisely, the 
algebraic quantum inverse scattering method permits us write the following equation 

t(X\z)$(X h ...X p ) = A(A, A 1; ■ ■ • , \ p \z)<S>(\ u ■ ■ ■ , X p ) - £ p^- (1-4) 

a=l A ~ A « 

Here t(X\z) denotes the transfer matrix of the rational vertex model in an inhomo- 
geneous lattice acting on an iV-fold tensor product of SU(2) representation spaces. 
$° = $(Ai, • • • X a -i, A, A Q+ i, A p ). JF q (Ai, • • • , Xp\z) and A(A, Ai, • • • , X p \z) are c num- 
bers. The vanishing of the so-called unwanted terms, T a = 0, is enforced in the usual 
procedure of the algebraic Bethe ansatz by choosing the parameters Ai, X p . In this 
case the wave vector $(Ai, • • • , X p ) becomes an eigenvector of the transfer matrix with 
eigenvalue A(A, Ai, ■ • • , X p \z). If we keep all unwanted terms, i.e. T a ^ 0, then the wave 
vector $ in general satisfies the equation (|1.4j) , named in || as off-shell Bethe ansatz 
equation (OSBAE). There is a neat relationship between the wave vector satisfying the 
OSBAE (|1.4j ) and the vector- valued solutions of the KZ equation (|1.3|) : The general 
vector valued solution of the KZ equation for an arbitrary simple Lie algebra was found 



by Schechtman and Varchenko [10]. It can be represented as a multiple contour integral 



. . . , z N ) = I ■ ■ ■ I X{X X , A p |z)0(Ai, X p \z)dX 1 ■■•dX p (1.5) 

The complex variables Zi, z^ of ( |1.5|) are related with the disorder parameters of the 
OSBAE . The vector valued function 0(Ai, X p \z) is the semiclassical limit of the wave 
vector $(Ai, X p \z). In fact, it is the Bethe wave vector for Gaudin magnets but 
off mass shell. The scalar function X(X\, ...,X p \z) is constructed from the semiclassical 
limit of the A(A = Zk] Ai, X p \z) and JF Q (Ai, ■ ■ • , X p \z). This representation of the Ap- 
point correlation function shows a deep connection between the inhomogeneous vertex 
models and the WZNW theory. 

In this paper we apply the Babujian-Flume ideas for osp(l\2) rational solution of 
the graded version of the YB equation 

The paper is organized as follows. In Section 2 we present the algebraic structure 
of the osp(l\2) vertex model. The inhomogeneous Bethe ansatz is read of from the ho- 
mogeneous case previously known j0|. We derive the Off-shell Bethe ansatz equations 
for the fundamental representation of the osp(l\2) algebra. In Section 3 , taking into 
account the semiclassical limit (77 — > 0) of the results presented in the Section 2, we 
describe the algebraic structure of the corresponding Gaudin model. Using the osp(l\2) 
Gaudin algebra, these results are extended for highest representations. In Section 4, 
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data of the Gaudin off-shell Bethe ansatz equations are used to construct solutions of 
the KZ equation. Conclusions are reserved for Section 5. 

2 Structure of the osp(l|2) Vertex Model 

We recall that the osp{l\2) algebra is the simplest superalgebra and it can be viewed 
as the graded version of sl(2). It has three even (bosonic) generators H, generating 
a Lie subalgebra sl(2) and two odd (fermionic) generators , whose non-vanishing 
commutation relations in the Cartan-Weyl basis reads as 

[H^X^] = ±X ± , [X + ,X~]=2H 

[H,V±\ = ±\v ± , [X ± ,V^] = V ± , [X ± ,V ± ] = 

{V*,^} = ±±X\ {V\V-} = - l -H (2.1) 
The quadratic Casimir operator is 

C 2 = H 2 + ^{X+, X-} + [V + , V-] (2.2) 

where {• , •} denotes the anticommutator and [• , •] the commutator. 

The irreducible finite-dimensional representations pj with the highest weight vector 
are parametrized by half-integer s = j/2 or by the integer j e N. Their dimension 
is dim(pj) = 2j + 1 and the corresponding value of C 2 is j(j + l)/4 = s(s + 1/2), 
5 = 0,1/2,1,3/2,... 

The representation corresponding to s = is the trivial one-dimensional represen- 
tation. The s > 1/2 representation contains two isospin multiplets which belong to 
isospin s and s — 1/2, denoted by \s,s,m) and \s,s — 1/2, m), respectively. The first 
quantum number characterizes the representation and the second and third quantum 
numbers give the isospin and its third component. After a convenient normalization 
of the states , a given s-representation is defined by 

H\s,s,m) = m\s,s,m) , 
X ± \s, s,m) = \J(s^fm)(s±m + 1) \s, s, m ± 1) 

V ± \s,s,m) = ±^(sTm) \s,s- 1/2, m± 1/2) 

H\s,s — 1/2, m) = m\s,s — 1/2, m) 
X ± \s,s- 1/2, m) = y/(s - 1/2 =f m)(s - 1/2 ± m + 1) \s, s - 1/2, mil) 

V ± \s,s-l/2,m) = ±^y/(s - 1/2 ± m + 1) \s, s, m ± 1/2) (2.3) 
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The fundamental representation has s = 1/2 and is given by 
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(2.4) 
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In (|24) the basis is 
considered as even and the second as odd, i.e., our grading is BFB. 
In the j-representation, the odd part has the form [12 



|, |, — I >. The first and third vectors will be 
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(2.5) 



where 

(V^V^V^,...^) 

(i>:,-i>;, :.u;, 2 u- , , i 

For the even part we can see from 
due to isospin integer: 



2 Wi> VT, \/j - 1, V2, VT, v'./ ■ 
1 



(- Jj, v 7 !, -Vj " 1, -VT, VJ')- (2-6) 



that H is diagonal and always has eigenvalue 



H = -diag(j, j - 1, 1,0, -1, -j). 



(2.7) 



Moreover, X ± are given by the si (2) composition which results in a clear relation with 
the odd part: X ± = ±4(y ± ) 2 . 



2.1 Graded Quantum Inverse Scattering Method 

Consider V — Vq © V\ a Z 2 -graded vector space where and 1 denote the even and odd 

s 

parts respectively. The components of a linear operator A ® B m. the graded tensor 

s 

product space V ® V result in matrix elements of the form 

(A & = (-)P(«(P(«)+PW) AaryBps (2.8) 
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and the action of the permutation operator V on the vector |a) 
given by 



e v <g> v is 



P I a) 



\p(a)p(/3) 



® | a) 



a/3 



{-) p{a)m S aP 5,6 (2.9) 



where p(a) = 1 (0) if \a) is an odd (even) element. 

Besides 1Z we have to consider matrices R = VIZ which satisfy 

Ri 2 (X)R 23 (X + n)R u (n) = R23^)Rn(X + /i)^ 23 (A). 



(2.10) 



The rational solution of the graded YB equation for the fundamental representation 
of osp(l\2) algebra was found by Kulish in |13|| . It has the form 
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(2.11) 



A + 77, 
A-77 + 

7] + 



x 2 (X,rj) = eiA 
Xr] 



x 3 {X,r]) = A - 



Xr) 



X + §77 



A + §77' 



At/ 



A + |V 



2s(A,T7) = 77, X 6 (A,77) 

ye(A, 77) = -x 6 (A,77), 



-e 2 



A77 



A + §7/ 



(2.12) 



where 

Xl(A,77) 
X 4 (A,7/) 

x 7 (X,r}) 

where e« = ±1, i = 1,2. Here we have assumed that the grading of threefold space 
is p(l) = p(3) = and p(2) = 1 and we will choose the solution of (|2.12|) for which 
ei = e 2 = 1. 

Let us consider the inhomogeneous vertex model, where to each vertex we associate 
two parameters: the global spectral parameter A and the disorder parameter z. In this 
case, the vertex weight matrix 1Z depends on A — z and consequently the monodromy 
matrix will be a function of the disorder parameters Zi. From now on we use a compact 
notation for the arguments with the shifted spectral parameter: (A|z) = (A — Z\, A — 
z N ). 
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The graded quantum inverse scattering method is characterized by the monodromy 
matrix T(X\z) satisfying the equation 



R{X - y) 



T{X\z)®T{^\z) 



T{^\z)®T{X\z) 



R(X-fi), 



(2.13) 



whose consistency is guaranteed by the graded version of the YB equation ( |2~ToD . 
T(X\z) is a matrix in the space V (the auxiliary space) whose matrix elements are 
operators on the states of the quantum system (the quantum space, which will also be 
the space V). The monodromy operator T(X\z) is defined as an ordered product of 
local operators C n (Lax operator), on all sites of the lattice: 

T{X\z) = C N {X\z)C N ^{X\z) ■ ■ ■ C x {X\z). (2.14) 

The Lax operator on the n th quantum space is given the graded permutation of (p.ll ): 
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(2.15) 



Note that L^(X), a, (3 = 1, 2, 3 are 3 by 3 matrices acting on the n th site of the lattice. 
It means that the monodromy matrix has the form 



where 



T i3 {X\z) 



hJ 



T(X\z) 



3 

E 

fcl,...,fcjV- 

1,2,3. 



( A x {X\z) B x {\\z) B 2 {X\z) 
C y {X\z) A 2 (X\z) B ? XX\z) 
\C 2 (X\z) C 3 (X\z) A 3 (X\z) 



(2.16) 



(2.17) 



The vector in the quantum space of the monodromy matrix T(X\z) that is annihi- 
lated by the operators Tij(X\z), i > j (Ci(X\z) operators, i = 1,2,3) and it is also an 
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eigenvector for the operators Ta{X\z) ( Ai(A|z) operators, i = 1, 2, 3) is called a highest 
vector of the monodromy matrix T(X\z). 

The transfer matrix t(A|z) of the corresponding integrable spin model is given by 
the supertrace of the monodromy matrix in the space V 

r(X\z) = ]T(-l) p(i) T u (X\z) = Ai(X\z) - A 2 (X\z) + A 3 (X\z). (2.18) 
i=i 

A detailed exposition of the graded quantum inverse scattering method can be 



found in the references [P^| ,p6| and |p7 | 



2.2 Inhomogeneous Bethe Ansatz 

Algebraic Bethe ansatz solution for the inhomogeneous osp(l\2) vertex model can be 
obtained from the homogeneous case. The only modification one have to do is local 
shifting of the spectral parameter A — > X — Z{. The algebraic Bethe ansatz solution 



of the homogeneous rational osp(l\2) vertex model was obtained by Martins in [jl8 
However, we will consider the rational limit of the trigonometric case presented in [19 
In the sequence we will work with some functions which will be now defined: 

Xi(A) _ A + 77 x 3 (X) _2X + r] 

Z \ A ) ~ — 77T — — x — j V\ A ) 



cu(X) 



x 2 (X) X y e (X) 2 V 

xi(A)a? 3 (A) _ 2A + 77 

x 4 (X)x 3 (X) - x 6 (X)y 6 (X) 2X-r] 



z{Xk Xj) ~ { l[xl - A MXj - X k ) if k< J j (2 - 19) 

We start defining the highest vector of the monodromy matrix T(X\z) in a lattice of N 
sites as the even (bosonic) completely unoccupied state 



|0> = ®£Li I (2.20) 




Using fl2.17| ) we can compute the action of the monodromy matrix entries on this state 



A(X\z)\0) = Xi{X\z)\0), Ci{X\z) |0) = 0, B t (X\z)\0) ^{0,|0)} 



N 



Xi(X\z) = l[x l (X-z a ), i = 1,2,3 (2.21) 



0=1 



Unlike to the simple Lie algebras cases (|1.4j) , the off-shell Bethe ansatz equations 
for the osp(l\2) algebra present in a more complicated form 
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n n j — 1 

r(A|z)*„(A!, A n ) = A*„(Ai, A n ) - £ jf>¥* B) + £ E *& n) *& ( 2 -22) 

j=l j=2 i=l 

where the Bethe vectors are defined as normal ordered states ^ n (Ai, • ■ ■ , A n ) which can 
be written with aid of a recurrence formula ]20 : 



* n (Ai, A n ) = Si(Ai)* n _i(A 2 , A n ) 

-B 2 (X 1 )J2 M l(Ail x\ ft ^(A fe -A i )^ n _ 2 (A 2 ,...,A i ,...,A n ) (2.23) 

i=2 ?/l A l _ A JV fc=2,fc# 

with the initial condition = |0) , l I / i(Ai) = £>i(Ai) |0). Xj denotes that the rapidity 

\j is absent: *(Ai, Xj, • • • , A n ) = ^(Ai, A^-i, • • • , X n ). 

Let us now describe each term which appear in the right hand side of the OSBAE 
for the osp(l\2) model ( for more details the reader can see fll9|): In the first term the 
Bethe vectors ( |2.23| ) are multiplied by c-numbers A = A(A, A 1; X n \z) given by 

" z{X-X k ) , vn| " x 2 (\-X k ) 



A = X^Xlz) J] 4h - A) - H n X 2 (X\z) J] " + X 8 (A|*) n 



fc=i^(A-A fe ) fc=1 s3(A-A fc ) 

(2.24) 

The second term is a sum of new vectors defined as 



*W = ( x2(A--A) gl(A|z) + ^V^T 5 3(Ak) ] ^n-i(Ai) (2.25) 

r(n) 



V.r.sA, A)"' ; ' y(A A ; ) 
and jFj"' 1 are scalar functions given by 

n n 

= X 1 {X j \z) J] Z{\ k - Xj) + (-) n X 2 (Xj\z) n Z{\ - A fc ). (2.26) 



Finally, the last term is a coupled sum of a third type of vector-valued functions 

, A A 

with intricate coefficients 



= B 2 (X\z)^ 2 (ii, Xj) (2.27) 



4 n) = G jl X 1 {X l \z)X 1 (X j \z) I] 2(A fc - A,)Z(A fc - Aj-) 

k=l,ky£j,l 

n 

-(-) n F, 7 X 1 (A / |z)X 2 (A,» n - - A fe ) 

k=l,k=/=j,l 
n 

-(-YFflX^X^X^z) I] Z[\ { - \ k )Z{\ k - Xj) 

k=l,kjtj,l 

n 

+H jl X 2 {X l \z)X 2 {X j \z) I] Z{Xj - X k )Z{X l - A fc ) (2.28) 
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where Gji , Yj\ , Fji and Hji are additional ratio functions defined by 
x 7 (\i -A) f z(Xi - A) x 5 (\j - A) 1 



H 



ji 



xs(Xi - A) y(Xi - Xj) lu(Xi - A) x 2 (Xj - A) y(X - Xi) 
y 7 (X - A/) f y 5 (A-A ; ) 1 

£ 3 (A - Xi) y(Xi - Xj) x 3 (A - A;) y(X - Xj) 

Y 1 { (X \) y ^ X ~ X ^ y 5 (X - Xi) y^ - Xj) 

'' y(A-A0l l l) x 2 (X-X J ) X2 (X-X l )x 2 (X l -X j ) 
Z/s(A - Aj f ysCAj-Aj) 1 + z(A - Aj 1 
x 2 {X - Xi) \ x 2 (Xi - Xj) y(X - Xt) u(X - X t ) y(X - Xj) 
Z/ 5 (A-A ; ) 1 \ 
x 2 {X - Xi) y{Xi - Xj) J 



F. 



(2.29) 



In the usual Bethe ansatz method, the next step consist in impose the vanishing of 
the so-called unwanted terms of ( |2.22j) in order to get an eigenvalue problem for the 



transfer matrix: 

We impose = , which also implies vanishing . Hence, \l/ n (Ai, A n ) is 
an eigenstate of t{X\z) with eigenvalue A , provided the rapidities Aj are solutions of 
the inhomogeneous Bethe ansatz equations 

" ( Xj-z a + V \ +1 " (A J -A fc + r / )(A fc -A J +r//2) 

L\\ Aj - z a ) 1 ) l^x^x.-^x.-x.-^m 

j = l,2,..,n (2.30) 



These equations were derived for the first time in W2 



3 Structure of the osp(l|2) Gaudin Model 

In this section we will consider the theory of the Gaudin model. To do this we need to 
calculate the semiclassical limit of the results presented in the previous section. 

In order to expand the matrix elements of T(X\z), up to an appropriate order in 77, 



we will start by expanding the Lax operator entries defined in (|2.15| ): 



r (?i)/\| \ -1 1 2if n 3 2^H^ -\- H n 3 
Ln{X\z) = I + 77- + 7^77 + 

A — Z n Z [A — Z n ) 



2 ' (A - Zn y 
' A - z n 2 (A - z n ) 2 



LS(X\z) = l-V^ + lv 2 T n H : +o(v 3 ), (3.1) 
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and for off-diagonal elements, we have 

L%\\\z) = ~V^ + o(v 2 ), L£\\\z)=r ] ^- + otf), 

A Z n A Z n 

Z&\\\z) = V^ + otf), L ( it(\\z)=r ] ^ + otf), 

-A — Z n A Z n 

L%\\\z) = v -^+o( V 2 ), L^(X\z)=r ] -^ + o(r ] 2 ). (3.2) 

A — Z n A — Z n 

Here H, and are matrices 3 by 3 given by fl2.4| ). 

Using ( |3.1| ) and ( |3.2| ) is now easy to derive the semiclassical expansion of the mon- 
odromy matrix entries: 



A 1 (X\z) 



( N T-t 

r(A|z) {i + vEt^- 



, „ 2 v | K®tt» + *-®*+-v-®y+ , 3 " [H 2 -H] a , , 3 x 

^ ^ (A-z )(A-z 6 ) + 4 a tl (A-z a ) 2 l + {V) 



a<b 

A 2 (X\z) = r(A|z){l + 7?0 



(A-z a )(A-z 6 ) 2 a ^(A-^l l/; 



a<6 



A 3 (A|z) 



77 ^ 1 (A - *»)(A - 4 ti (A - ' ^ J 



a<6 



f ^ v- 

^i(Ak) = -B 3 (X\z) = T(X\z)\- v Y,t^ + o(v 2 ) 



a=l A - Za 



Ci(A|*) 



c s (A|z) = r(A|z) Lf;J2_ + (,ft] 

I a=l A - Z a J 



iv v- 

A„ , 2 



S 2 (A|z) = r(A|*)^£ r f- + fo 



a=l 



C 2 (A|*) = r(A|s) f7 5]_2_ + (^) (3.3) 

I a=l A ~ Za ) 

where W = 2H, X ± = 2X ± , V ± = 2V ± and T(\\z) = Ua=i( x - z a) is a common 
factor which can be absorbed after a convenient normalization. From the definition 
Q2.1S| ) it follows that the semiclassical expansion of the normalized transfer matrix 
t(\\z) = t(X\z)/T(X\z) has the form 

[l]a 



{N (~t q N 



+ o(rf ) (3.4) 
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where 



N 



G ^ Hg®H h + \{X+k>Xf + X;£> X b + ) + V+ I V£ - Yl ® yl 

b^a Za ~ Zb 

Here we have used the symmetry Gkj = VjkGjkVjk = Gjk and the identity 



1 



1 



1 



(3.5) 



(3.6) 



(X-z a )(X-z b ) \X-z a X-z b Jz a -z b 

The Gaudin model is defined by the residue of t(X\z) (|3.4j ) at the point A = z a . 
This results in N non-local magnets whose Hamiltonians G a , a = 1, 2, N are given 
by (p.5|) and satisfied the property J2 a =i G a = 0. 



3.1 osp(l|2) Gaudin Algebra 

Let us now consider the semiclassical limit of the fundamental commutation relation 
( [2.13| ). Let us first write the i?-matrix (|2 . 1 1 ) in the form [13|: 



R(X, V )=r]l+XV + 



Xrj 



X + \r) 



U 



(3.7) 



where X is the identity operator, V is the graded permutation operator (|2.9|) and U 
is the rank-one projector U 2 =U. Using the normalization of the previous section we 
can write the semiclassical expansions of T and R in the following form 



T{X\z) 



l + r}L(X\z) + o(r] 2 



R(X-fj) = V 1 + r]r(X — fJL) + o(r] 2 ) 



From 



we can see that the "classical L-operator " has the form 
L{X\z) = 



H(X\z) -V-(X\z) X~(X\z) 
V + (X\z) V-(AI-z) 
X+{X\z) V+{X\z) -H(X\z) 



(3.8) 



(3.9) 



where 



N 



nj N ■> l- 

Efj. v-(A|,) = i;Ji-, 

a=l A - Z a a=l A ~ Z a 

0=1 A _ Z a a =l A ~ Z a 



H(X\z) 
V + (X\z) 



N 



x-(x\ Z ) = y: 



x: 



o=l A - z a 



(3.10) 



and the classical r- matrix is obtained by the semiclassical expansion of (3.7). It has 
the form 



r(A- n) 



1 

A — n 



(V-U) 



(3.11) 
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Substituting ( |3.9| ) and ( |3. 10| ) into the fundamental relation ( [2.13] ) we will get, from 
the first non trivial identity , the following equation: 



VL{X\z) <g> L(fj,\z)+Vr{\-ti) 



L(X\z) <g> 1 + 1 ® L(fi\z) 



L{n\z) ® L{\\z)V 



L(\\z) ® 1 + 1 <g> L(fj,\z) 



Vr(\-n) (3.12) 



whose consistence is guaranteed by the graded classical YB equation. 

From ( |3.12| ) we can derive commutations relations between the matrix elements of 
L(X\z). They are the defining relations of the osp(l\2) Gaudin algebra presented in 

0: ' 



[H(\\z),H(n 
[V + (X\z),X + (fi 

[H(X\z),V-(fi 
[H(\\z),X-(n 
[H{\\z),V + {n 
[H{\\z),X + (im 
[X-(X\z),X + (ji 



\y-{\\z),x-{»\z)] = [X-(\\z),X-(ti\z)\ = 
[X + (\\z),X + (ji\z)]=0 



-2[X-(\\z),V + (fi\z)\ 

2{V-(A|z),V-(Ai|s)} = 
-2[V-(X\z),X + (i,\z)] 
2{V + (X\z),V + ^\z)} = 
4{V-(A|z),V + (/i|z)} = 



A — fj, 
2 



[V-(\\z)-V-(»\z)} 



A — jj, 



1 



A — n 
2 



[X-(X\z)-X-(fi\z)\ 
[V + (X\z)-V + (»\z)] 



[X+(\\z)-X+(v\z)} 



A — [A 

' [H(X\z)-H(n\z)] (3.13) 



A — jj, 



It follows from ( |3.12| ) that the Gaudin algebra is the semiclassical limit of the com- 
mutation relations of the matrix elements of the monodromy matrix used in the usual 
Bethe ansatz method. 

A direct consequence of these relations is the commutativity of t(X\z) 



[t{X\z),t{v\z)} = 0, VA,/i 



(3.14) 



and as consequence it follows the commutativity of the Gaudin Hamiltonians G a ,a = 
l,..,iV. 

3.2 Gaudin Off-shell Bethe Ansatz Equations 

In order to get semiclassical limit of the OSBAE fl2.22|) we first consider the semiclassical 
expansions of the Bethe vectors defined in ( |2.23| ), ( |2.25| ) and ([2.27]) : 



*n(A ls ...,A n ) = (-r ] ) n ^ n (X 1 ,...,X n \z) + o(r ] n+1 ) 
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4 J 

= -(-r ? ) n - 1 A'-(A|z)$„_ 2 (A 1 ,...,A i ,A„...,A ri ^) + o(r/ n ) (3.15) 



where 

$„(A 1 ,...,A n |z) = V~(Ai|z)$„_i(A 2) ...,A n |z) 

-X-{\ x \z) J2 T^v$n- 2 (A 2 , Aj, K\z) (3.16) 
j=2 Al ~ A i 

with $ = |0) and $i(Ai|z) = V~(Ai|z) |0). 

Here we would like make a few comments on the structure of these vector-valued 
functions. In (|3.16|) they are written in a normal ordered form. Since we are working 
with fermionic degree of freedom, the function $ n (Ai, ...,X n \z) is totally antisymmetric. 

^n(Ai, Aj_i, Aj+i, Aj, X n \z) = — $ n (Ai, Aj_i, Xi, Aj+i, X n \z) (3-17) 

To see this one can use the Gaudin algebra ( |3.13| ). For instance, in its antisymmetric 
form the Bethe vector $ 2 reads as 

$ 2 (Ai,A 2 |*) = ^[V-(A 1 |^)V-(A 2 |^-V-(A 2 |z)V-(A 1 |z)]|0) 

_i [ ^±^m m (3 , 8) 

Next we consider the semiclassical expansions of the c- number functions ( |2.24 ), 
( g^g ) and (gg|) presents in the OSBAE ( ^22|) 

A = 1 + 2r? 2 A (2) + o{rf) 
Jf> = Vi-yf^ + oiv 2 ), 4 n) = 2r ? 3 (-)^F^ + o( ?7 4 ) (3.19) 



where 



and 



N i / N 1 \ n I / n i 

a (2) = F — - — F + v — - — y 

tl ( A - *0 ^ * - *>) tl ( A " Afc ) Vfe Afc - A i , 
N n -1 i q AT -i 

+ EE A z XT a + I^o^ (3 ' 20) 

o=l fc=l A ~ Za Afc ~ A 4 a=l l A ~ 



/i n) = eA-e 



A i - A fc a =i A i ~~ ^ 

(n) * An) 



p{n) _ _Jj_ l _ Jl ; 21) 

^ A — Xj Xj — Xi X — Xi A/ — Aj 

13 



Substituting these expressions into the ( |2.22| ) and comparing the coefficients of the 
terms 2(— r]) n+2 we get the first non-trivial consequence for the semiclassical limit of 
the OSBAE: 



E T^^(Ai, A n ) = .V i: <I>„(A,. .... A n ) + E ' " 



where 



f\ X ~ z a 

&n = V-(\\z)<$> n „ 1 (\ 1 ,...,\ j ,...,\ n \z) 



i= x A -A, 



(3.22) 



-x-(\\z)y: 



A A 

\i, \j, X n \z) 



(3.23) 



with I— I + 1 for I < j or /= I for / > j. 

Here we observe that the OSBAE fl3.22j ) is similar to the OSBAE presented by 
Babujian and Flume for simple Lie algebras. This result could be expected since the 
superalgebra osp(l|2) has many features which make it very close to the Lie algebra 



Now, we take the residue of ( |3.22j ) at the point A = z a to get the Gaudin off-shell 
equations: 

Ji 



G $„(Ai, A n |z) = # a $ n (Ai, A„|z) + E 



1=1 Z a ~ \ 



1,2, ...,N 



where /j is given by (|3.21|) , g a is the residue of A^ 2 ) ( |3.20|) at A 



N 



b^a Za Zb 1=1 Za Al 



E 



(3.24) 



(3.25) 



A A 
$>n-2{\l, Xj \z) 



(3.26) 



and 4> l n is the residue of <3>^ ( E>.23| ) at A = z a 

ft = VrQn-t&j \z) -xfjr -t^- 

In this way we arrive to one of the main problem solved by the theory of the Gaudin 
model, i.e., the determination of the eigenvalues and eigenvectors of the commuting 
Hamiltonians G a ( |3.5|) : g a is the eigenvalue of G a with eigenfunction $ n provided A; 
are solutions of the Bethe ansatz equations 

1 A 1 



fi n) = o =► E 



N 

E 



k^l X l ~ X k ,,=1 A; - Z a 



I = 1, 2, n. 



(3.27) 



These Bethe ansatz equations are the semiclassical limit of the inhomogeneous Bethe 
ansatz equations ( 2.3C ). 
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4 Knizhnik-Zamolodchickov Equation 



Since the spectra of the osp(l\2) Gaudin Hamiltonians are now known, the next step 
is to calculate correlation functions. 

Before we present the result of this section, let us observe a very important conse- 
quence due to semiclassical approach. While the solutions of the quantum YB equation 
( |1.1| ) depend essentially on the representation, the corresponding classical solutions 
fij(X) G gi <8> gj , g being a Lie algebra, can be written in an invariant form, i.e. inde- 
pendent of the representation. Therefore, one may extend the previous semiclassical 
results beyond the fundamental representation. To do this we first define a r-matrix 
and the corresponding L-operator for the next representation (s = 1) of the osp(l\2) 
algebra presented in (|2.3|). The classical r-matrix is constructed out of the quadratic 
Casimir in a standard way |T5[, 



r(A-/x) 



A — jJL 



[H 



®H+^(X" 



<§) X~ + X~ <§) x~ 



+ v+ ® v~ - v- ® v + ] (4.1; 



Using the basis ( j2.3|) with the grading BFBFB one can verify that this classical r-matrix 
satisfies the graded version of the classical YB equation (OF). To write (|3.11| ) in the 



form ( |4.1| ), we recall the symmetries of the YB solution ( |2.12| ) by considering the case 
€i = 1 and e 2 = — 1. 

The classical L-operator is obtained from the semiclassical limit of the canonical 
identification of the Lax operator L with 7£-matrix [L4]]. In the s = 1 representation 
it has the form 



L(X\z) 



V 



2H(\\z) 
-V2V + (X\z) 
V2X + (X\z) 







V2V~(X\z) 
H(\\z) 
V + (X\z) 
X+(X\z) 




y/2X-(X\z) 
V-(X\z) 


-V + (A|.) 
V2X + (X\z) 





X-(X\z) 
V-(X\z) 
-H(X\z) 
V2V+(X\z) 








V2X~(X\z 
V2V-{\\z 
-2H(X\z 



) / 



(4.2) 

Substituting ( |4.2j ) and ( fCT| ) into the fundamental relation ( B.12| ) we will get the 
same defining relations of the osp(l\2) Gaudin algebra ( pj,13| ). 

The fundamental property of the operators ( p.l(J[ ) is that they form the highest 
weight module of the infinite-dimensional Lie superalgebra. As in the s/(2) case pre- 
sented by Sklyanin in |E3 , it is characterized by the vacuum |0) , its dual (0| 



H(X\z) |0> 
(0\H(X\z) 



h(X\z) |0), 
h(X\z) (0|, 



V + (X\z)\0) 
(0\V-(X\z) 



0. 
0. 



X+(X\z)\0) 
(0\X-(X\z) 








(4.3) 
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and the highest weight functions h(\\z) 



N 2s , 

#) = ErT- Sa = -,i,-,... (4.4) 



In this way we are going beyond the fundamental representation s a = |. 

The semiclassical expansion of t(X\z) in the s-representation is easily obtained from 

(B 

t(X\z) = 1 + 2, 2 {£ + £ + o(, 3 ) (4.5) 

U=l A - a= i (A - 2 a J J 

Here we remark that this result is general. It comes from the structure of the semiclas- 
sical method, where t(X\z) is defined as supertrace of \L 2 (X\z) and the second sum in 
Q4.5| ) is identified with the Casimir of the subalgebra si (2). Consequently, the Gaudin 
off-shell equations in the s-representation are known. They are again given by ( |3.24| ), 
where the coefficients and g a are 

kf L t M - A fc a=1 A t - z a 

N n 

9a = 4£^-2£^Y (4.6) 

b^a Z a Z b i=i z a *l 

and $ n (Ai, A n ) for s-representations again has the form as Eq.( p.l6|) . 

Let us now define the vector-valued function zn) through multiple contour 

integrals of the Bethe vectors ( p.!6|) 

...,z N )=f---f X(X\z)^ n (X\z)dX 1 ...dX n (4.7) 

where X (\\z) = X (Ai, A n , z±, z^) is a scalar function which in this stage is still 
undefined. 

In analogy with the sl{2) case M, we assume that *&(zi, z^) is a solution of the 
equations 

=G a *( Zl ,...,z N ), a = 1,2,..., iV (4.8) 

oz a 

where G n are the Gaudin Hamiltonians 03.51) and k is a constant. 



By construction, these equations are the graded version of the KZ equations f25| , 
26| . Hence, we can interpret ^(zx, ■■■,Zn) as a integral representation for the iV-point 



correlation function in an chiral conformal field theory with osp(l\2) symmetry [27], |28 
Substituting (O) into (ED we nave 



K — ^ = J ' "/^ $w(A| ^ + ^ 1 ' }dXi...dX n (4.9) 
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From ( |3.16| ) and using the Gaudin algebra ( |3.9| ) one can derive the following identity 



<9$ r 



d 



i d\i \Xi — z at 
which allows us write (14.91) in the form 



(4.10) 



K 



dz a 



■ dxm (xlz) A, )lr dx(\\z) 



i=l 

W ^ (X{\\z)<fc 



, , <9A/ \ Xi - Z a 



]dXi...dX n 



\]dXi...dX n 

(4.11) 



It is evident that the last term of (|4.11|) is vanishes, because the contours are closed. 
Moreover, if the scalar function ^(A^) satisfies the following differential equations 



dX(X\z) 
« g z =g a X(X\z), 



(4.12) 



we are recovering the off-shell Bethe ansatz equations for the Gaudin magnets (|3.24 ) 
from the first term in ( [4.1 1| ) . Taking into account ( |4.6|) we can see that the consistency 



condition of the system ( |4.12| ) is insured by the zero curvature conditions dfj / d 
dg a /dXj. 

The solution of ( [4.1 2[ ) is 



N n N,n 

x(x\ z ) = ]j(z a - Zb y s ^i* H(Xj - A,) 1 /" H(z a - Xj)~ 2sa/K . 

a<b j<k a,j 



(4.13) 



This function determines the monodromy of ^(z\, z^) as solution of the osp(l\2) KZ 
equation ( |4.8| ) and these results are in agreement with the Schechtman-Varchenko con- 
struction for multiple contour integral as solutions of the KZ equation in an arbitrary 
simple Lie algebra [HI. 



5 Conclusion 

In the first part of this paper we have considered an integrable representation of the 
osp(l\2) algebra to study the algebraic Bethe ansatz for an inhomogeneous vertex 
model. 

In the second part, through the semiclassical limit of the quantum inverse scattering 
equations we got arrive to the theory of the Gaudin model. Here, with aid of the Gaudin 
algebra, we extended the results presented in the first part for highest representations 
and find the spectra of N commutating Hamiltonians. 
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We believe that the procedure presented here will work for all models for which 
the usual algebraic Bethe ansatz is known. Therefore, it should be interesting extend 
this procedure for all algebraic versions of the Bethe ansatz. In particular, the natural 
candidate will be the nested algebraic Bethe ansatz |29|. In ||30|| , a system of matrix 



difference equations was solved by means of a nested version of a off-shell Bethe ansatz. 

Regarding our results about the Appoint correlators in the osp(l\2) conformal field 
theory, there are several issues left for future works. For instance, we would like to con- 



nect this result with those presented in the osp(l\2) current algebra [31, 32]. It should 



be interesting since it appears in the quantization of two-dimensional supergravity in 



the light-cone gauge [Q. Moreover, the affine version of the osp(l\2) algebra is the 
starting point in the construction of the N = 1 superconformal minimal models by 
means of the Hamiltonian reduction procedure |3l 



Acknowledgment: This work was supported in part by Fundagao de Amparo 
a Pesquisa do Estado de Sao Paulo-FAPESP-Brasil and by Conselho Nacional de 
Desenvolvimento-CNPq-Brasil. 



References 

[1] R.J. Baxter, Exactly solved models in statistical mechanics, Academic Press-1982, 
London. 

[2] AA. Belavin and V.G. Drinfeld, Func. Anal. Appl. 16 (1982) 159. 

[3] V.G. Knizhnik and A.B. Zamolodchikov, Nucl. Phys. B247(1984) 83. 

[4] L.D. Faddeev and LA. Takhtajan, Uspekhi Mat. Nauk 34 (1979) 13. 

[5] V.E. Korepin , N.M. Bogoliubov and A.G. Izergin, Quantum Inverse Scattering 
Method and Correlation Functions, Cambridge-1993. 

[6] V.E. Korepin, Commun. Math. Phys. 94 (1982) 67-113. 

[7] H.M. Babujian and R. Flume, Mod. Phys. Lett. A9 (1994) 2029. 

[8] H.M. Babujian, J. Phys. A: Math. Gen 26 (1993) 6981. 

[9] B. Feigin , E. Frenkel and N.Yu Reshetikhin, Commun. Math. Phys. 166 (1994) 
27. 

[10] V.V. Schechtman and A.N. Varchenko, Invent. Math. 106 (1991) 139. 



18 



[11] M. Gaudin, J. de Physique Z7 (1976) 1086. 

[12] P.P. Kulish, Quantm Superalgebra osp(2\l) LOMI Preprint-1988. 

[13] P.P. Kulish, Lett. Math. Phys. 10 (1985) 85. 



[14] L.D. Faddeev, How Algebraic Bethe Ansatz works for integrable model, [hep- 
| th/ 9605187 . 



[15] P.P. Kulish and K. Sklyanin, J. Soviet Math. 19 (1982) 1596. 
[16] P.P. Kulish, J. Soviet Math. 35 (1985) 2648. 
[17] F.H.L. Essler and V.E. Korepin, Phys. Rev. B46 (1992) 9147. 
[18] M.J. Martins, Nucl. Phys. B450 [FS] (1995) 768. 

[19] A. Lima-Santos, Bethe Ansatze for Nineteen Vertex Models, (iep-th/9807219 
[20] V.O. Tarasov, Teor.Math.Phys.76 (1988) 793. 
[21] T. Brzezinski and A.J. Macfarlane, J. Math. Phys. 35 (1994) 3261. 
[22] M.J. Martins and PB. Ramos, Nucl. Phys. B500 [PM] (1997) 579. 
[23] L. Frappat , P. Sorba and A. Sciarrino, Dictionary on Lie Superalgebras, [hep- 



th/9607161| . 



[24] K. Sklyanin, Generating function of correlators in the si? Gaudin model, |hcp- 
| th/ 9708007| . 

[25] P. Di Vecchia P, V.G. Knizhnik, J.L. Petersen and P. Rossi, Nucl. Phys. B253 
(1985) 701. 

[26] J. Fuchs, Nucl. Phys. B286 (1987) 455. 
[27] S. Nam, Phys. Lett. B187 (1987) 340. 

[28] E.B. Kiritsis and G. Siopsis, Phys. Lett. B184 (1987) 353 (Phys. Lett. B189 (1987) 
489). 

[29] H.J. de Vega, Int. J. Mod. Phys. A4 (1987) 2371. 

[30] H.M. Babujian, M. Karowski and A. Zepletal, J. Phys. A: Math.Gen. 30(1997) 
6425. 

19 



[31] L.P. Ennes, A.V. Ramallo A V and J.M. Sanches de Santos, Nucl. Phys. B491 
(1997)574. 

[32] L.P. Ennes and A.V. Ramallo, Nucl. Phys. B502 (1997) 574. 

[33] A.M. Polyakov and A.B. Zamolodchikov, Mod. Phys. Lett. A3 (1988) 1213. 

[34] M. Bershadsky and H. Ooguri, Phys. Lett. B229 (1989) 374. 



20 



